First-principles study of Ti-doped sodium alanate surfaces 
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We have performed first-principles calculations of thick slabs of Ti-doped sodium alanate 
(NaAlH4), which allows to study the system energetics as the dopant progresses from the sur- 
face to the bulk. Our calculations predict that Ti stays on the surface, substitutes for Na, and 
attracts a large number of H atoms to its vicinity. Molecular dynamics simulations suggest that the 
most likely product of the Ti-doping is the formation of H-rich TiAU (n > 1) compounds on the 
surface, and hint at the mechanism by which Ti enhances the reaction kinetics of NaAlH4. 
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One key issue for the advancement of fuel-cell technolo- 
gies is the development of safe and inexpensive ways of 
storing hydrogen aboard vehicles. In recent years sodium 
alanate (NaAlH4) has become one of the most promis- 
ing systems to achieve this goal. The kinetics of the re- 
versible reaction by which pure NaAlH4 releases hydro- 
gen is relatively slow. However, it was discovered that 
a few percent of Ti doping increases the reaction rates 
dramatically 1], bringing the system close to what is re- 
quired for practical applications. In spite of extensive 
investigations, the mechanism by which Ti enhances the 
reaction kinetics is still unknown In fact, even 

the location of the dopants, on the surface ^ or in the 
bulk of the system 0, 13 , and the reactions they cause re- 
main unclear. Some recent experimental results suggest 
the formation of a TiAla compound on the surface of the 
system 7J, but a consensus on this issue is yet to emerge. 

In this Letter we present a theoretical study address- 
ing the following outstanding questions. Where is Ti 
in doped NaAlIl4, on the surface or in the bulk? If Ti 
stays on the surface, does it substitute for Na or Al, and 
what are the reactions and structural rearrangements it 
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The calculations were performed within the general- 
ized gradient approximation to density functional the- 
ory as implemented in the code SIESTA 0. We used a 
localized basis set including double-^ and polarization or- 
bitals, and Troullier-Martins pseudopotentials We 
tested the convergence of our calculations with respect to 
the k-point and real space meshes. For example, for the 
001-surface calculations involving 4500 A^ supercells and 
about 162 atoms, we used a 2x2x1 k-point grid and a 
150 Ry cutoff for the real space mesh. In our SIESTA cal- 
culations some accuracy is sacrificed, mainly because of 
the small basis set we use, so that bigger systems can be 
studied. Hence, we checked some representative SIESTA 
results with a more accurate ab initio method, using a 
plane-wave basis set and ultrasoft pseudopotentials [ll| 
as implemented in the code CASTEP 1^. 

We started by determining which doping models of 
NaAlH4 surfaces are most energetically favorable. For 
that purpose we considered the two slab-type supercells 
in Figs, ^-b and QJ-d, which correspond to the most 
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(b) 1 00-slab (side view) 
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FIG. 1: Top and side views of the 100-slab (a-b) and the 001- 
slab (c-d) supercells of NaAlH4. Small and large gray spheres 
represent H and Na atoms, respectively. The Al atoms are at 
the center of the tetrahedra. 



natural terminations of tetragonal NaAlH4 [along the 
(100) and (001) directions, respectively]. Note that these 
supercells can be respectively regarded as composed of 
3x3x4 and 3x3x3 NaAlH4 groups (i.e., 216 and 
162 atoms). In both cases we imposed that the atoms in 
the deepest layer be fixed at the bulk atomic positions, 
and allowed about 13 A of empty space between slabs. 
The bulk structural parameters (a = 6 = 5.01A and 
c= 11.12 A) were obtained from first-principles and are 
in good agreement with previously published results Jl^ • 
We considered a structure to be relaxed when residual 
force components are smaller than 0.03 eV/A. NaAlH4 
is a fairly ionic system T^, TE^ and the pure surfaces do 
not present particularly relevant features. The electronic 
energy minimizations converge easily and the structural 
deviations from the bulk are relatively small. 

In order to investigate the Ti doping of these surfaces, 
we computed the cohesive energy (£^coh) of a number of 
experimentally motivated doping models the sub- 
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TABLE I: Calculated cohesive energies, given in eV and per 
NaAlH4 formula unit, for several doping models (see text). 
The column marked with an asterisk shows CASTEP results. 
In all the cases, ^coh of the pure surface is taken as the zero 
of energy. The pure ^coh values are given in the last line. 
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Doping model 


(100) 


(001) 


(001)* 


Ti @ Al 


0.05 


0.02 


0.04 


Ti @ Al + Na^ 


-0.01 


-0.09 


-0.06 


Ti @ Al + 2Na^ 


-0.11 


-0.28 


-0.25 


Ti @ Na 


0.11 


0.12 


0.13 


Ti @ Na + Na^ 


0.04 


0.03 


0.04 


Ti @ Na + 2Na^ 


-0.01 


-0.07 


-0.05 


undoped 
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20.15 
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stitution of Al and Na by Ti, as well as the occurrence of 
Na vacancies near the dopant. We use, for example, the 
notation "Ti@Na+Na^" to refer to the case in which Ti 
substitutes for Na and there is a Na vacancy next to it. 
In all the calculations we considered only one Ti dopant 
in the supercell. 

The cohesive energies shown in Tabled were obtained 
as the sum of the individual atom energies minus the 
energy of the system. For both surfaces, £^coh of the pure 
case is taken as the zero of energy. Hence, the results 
in the table give minus the energy change in reactions of 
the form Ti+Na27Al27Hio8 ^ Al+Na27Al26TiHio8, and 
measure the relative stability of pure and doped systems. 
A positive (and large) entry in the table indicates that, 
in principle, it is feasible to obtain the doped surface. 

We find that both Ti@Na and Ti@Al are energetically 
more stable than the pure surfaces, i.e., the system gains 
energy by accepting a Ti dopant and releasing a Na or Al 
atom. In addition, Ti@Na is found to be the most favor- 
able substitution, and the differences between the (100) 
and (001) surfaces are minor. Interestingly, the most fa- 
vorable surface doping models turn out to be the same as 
in the bulk case [l^ . We also find that the energy cost of 
forming Na vacancies is relatively small, suggesting they 
are likely to occur in the real doped system. Finally, Ta- 
ble ffl shows that the SIESTA and CASTEP results are in 
good agreement. 

The obtained structural relaxations closely resemble 
those described in Ref. 13] for Ti-doped bulk NaAlH4. 
In the cases where Ti substitutes for Na, the dopant drags 
neighboring H atoms towards itself and greatly stretches 
the corresponding Al-H bonds. The Ti for Al substitu- 
tion leads to an expansion of the TiH4 group, as Ti is 
larger than Al. This size difference is probably the rea- 
son that Ti@Na has a larger cohesive energy than Ti@Al, 
even if the nominal valences of the atoms would suggest 
otherwise. 

Next we studied whether the Ti dopant stays on the 
surface or, rather, penetrates into the bulk of the system. 
For that purpose, we considered slabs that are very thick 



FIG. 2: Calculated cohesive energies of various doping models 
(see text) as a function of the Ti dopant location (see right 
panel). In all cases the pure surface Ecoh is taken as the zero 
of energy. 



along the direction normal to the surface (see Fig.EJ, and 
computed the cohesive energy of the system as a function 
of the dopant position. More precisely, we worked with 
supercells formed by 7x2x2 and 2x2x7 NaAlH4 groups, 
respectively, for the (100) and (001) surfaces. The atoms 
in the deepest layer are fixed at the bulk positions. Such 
supercells allow us to study the dopants up to three layers 
away from the surface, which we found is enough to reach 
the bulk limit. We focused on the doping models found 
to be more stable than the pure surface, i.e., Ti@Na, 
Ti@Al, and Ti@Na-hNa^. 

Figure [2 shows the results. In all cases £^coh of the 
pure system is subtracted so that the value at the sur- 
face approximately coincides with that of Table ^ (The 
agreement between Tableland Fig. [21 should not be per- 
fect, since the respective systems differ in dopant con- 
centration.) Since a larger ^coh implies greater stability, 
we find that in all the cases considered it is energetically 
favorable for the Ti atom to remain on the surface of the 
system. In fact, the dopant always finds it favorable to 
locate in the outermost layer, except in the case of Ti@Al 
on the (001) surface. 

As mentioned above, we find pure NaAlH4 surfaces 
to be electronically and structurally similar to the bulk 
of the material. Hence, we think the dopants stay on 
the surface simply because this allows larger structural 
relaxations to accommodate them. 

Finally, in an effort to understand better the reactions 
and structural relaxations that occur in real NaAlH4 sur- 
faces, we performed ab initio molecular dynamics (MD) 
simulations at 300 K. We worked with the (001) slab de- 
picted in Fig.n^-d and studied the doping models Ti@Na, 
Ti@Na+Na^, and Ti@Na+2Na^. Our MD runs were 
about 2 ps long, which is hardly enough to thermalize 
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FIG. 3: Pair distribution function obtained from a 300 K ab 
initio MD simulation of the Ti@Na+2Na'' 001-surface (see 
text). 



a system as reactive as this. Nevertheless, we were able 
to study the average bonding of Ti by calculating pair 
distribution functions (PDF's; Fig. [^1 shows the result for 
Ti@Na+2Na^), and reached two distinct conclusions. 

In all the cases considered Ti is surrounded by a large 
number of H atoms. More precisely, the first Ti-H PDF 
peak is centered at about 1.9 A and corresponds to 7 hy- 
drogens for Ti@Na and 8 for Ti@Na+2Na^, Ti@Na+Na^ 
being an intermediate case. These hydrogens travel con- 
tinuously and freely between the dopant and the Al atoms 
they were originally bound to. The presence of the 
dopant thus results in a large number of loosely bound 
hydrogens that are relatively close to each other. This 
might well facilitate H2 formation and release, conse- 
quently enhancing the kinetics of the reactions in the 



decomposition/formation of doped alanates. Note that 
in Ref. 13] bulk Ti dopants were found to bind strongly 
a few H atoms at K, and that is probably also the case 
for surface dopants. However, the dynamical situation at 
300 K seems to be one in which many hydrogen atoms 
approach the dopant and none of them is actually tightly 
bound. 

The other effect that we observed in all our MD runs is 
the formation of Ti-Al bonds. The first Ti-Al PDF peak 
is centered at about 2.75 A and involves 2 Al atoms in the 
Ti@Na and Ti@Na+2Na^ cases, and almost 4 in the case 
of Ti@Na+Na^. Taking into account that the typical 
Ti-Al distance in TiAla systems is about 2.80 A |T^ . 
this clearly hints at the formation of some kind of TiAl^ 
(n > 1) compound on the surface of the real material. 
The calculations suggest such a compound should be rich 
in interstitial hydrogen. Unfortunately, our supercells 
are not big enough to study the size and structure of 
these hypothetical TiAl^ compounds. Nevertheless, our 
results clearly support the conclusions of the recent X-ray 
absorption study of Ref. gi] . 

In summary, we have carried out a first-principles 
study of Ti-doped NaAlH4 surfaces. Our conclusions are: 
(1) It is energetically favorable for the Ti dopants to stay 
on the surface rather than to penetrate into the bulk of 
the material. (2) Ti substitutes preferentially for Na, and 
formation of neighboring Na vacancies is likely. (3) The 
Ti dopant attracts to its neighborhood a large number of 
H atoms (up to 8). These hydrogens are loosely bound 
and close to each other, which might facilitate H2 for- 
mation and release. (4) We find strong indications that 
a TiAl^ (n > 1) compound forms on the surface of the 
alanate, in support of recent experimental reports 7]. 

We thank T.J. Udovic and C. Brown for their com- 
ments on the manuscript. 



[1] B. Bogdanovic and M. Schwickardi, J. Alloys Compd. 
253, 1 (1997). 

[2] K.J. Gross, G.J. Thomas, and C.M. Jensen, J. Alloys 

Compd. 330-332, 683 (2002). 
[3] T. Kiyobayashi, S.S. Srinivasan, D. Sun, and C.M. 

Jensen, J. Phys. Chem. A 107, 7671 (2003). 
[4] K.J. Gross, E.H. Majzoub, and S.W. Spangler, J. Alloys 

Compd. 356, 423 (2003). 
[5] G.J. Thomas, K.J. Gross, N. Y. C. Yang, C. M. Jensen, 

J. Alloys Compd. 330-332, 702 (2002). 
[6] D. Sun, T. Kiyobayashi, H.T. Takeshita, N. Kuriyama, 

and C.M. Jensen, J. Alloys Compd. 337, L8 (2002). 
[7] J. Graetz, J.J. Reilly, J. Johnson, A.Y. Ignatov, and T.A. 

Tyson, Appl. Phys. Lett. 85, 500 (2004). 
[8] J.P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. 

Lett. 77, 3865 (1996). 



[9] J.M. Soler et al, J. Phys.: Cond. Matt. 14, 2745 (2002). 
[10] N. Troullier and J.L. Martins, Phys. Rev. B 40, 2980 
(1993). 

[11] D. Vanderbilt, Phys. Rev. B 41, 7892 (1990). 

[12] M.C. Payne, M.P. Teter, B.C. Allan, T.A. Arias, and 

J.D. Joannopoulos, Rev. Mod. Phys. 64, 1045 (1992). 
[13] J. Iniguez, T. Yildirim, T.J. Udovic, M. Sulic, and C.M. 

Jensen, Phys. Rev. B 70, 060101 (2004). 
[14] A. Aguayo and D.J. Singh, Phys. Rev. B 69, 155103 

(2004). 

[15] A. Peles, J.A. Alford, Z. Ma, L. Yang, and M.Y. Chou, 

Phys. Rev. B 70, 165105 (2004). 
[16] C. Colinet and A. Pasturel, J. Phys.: Cond. Matt. 14, 

6713 (2002). 



